Numerical solution of buoyancy-driven flow problems in two spatial dimensions is presented. A high-order spectral method is applied for the spatial discretization, while the temporal discretization is done by operator ...